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ABSTRACT 

Context. The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently 
single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. 

Aims. One member of this group, the late O-type star VETS 399, is found to be unexpectedly X-ray bright for its bolometric luminos¬ 
ity - in this study we aim to determine its physical nature and the cause of this behaviour. 

Methods. To accomplish this we performed a time-resolved analysis of optical, infrared and X-ray observations. 

Results. We found VETS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demon¬ 
strates complex emission profiles in the lower Balmer series and select He i lines - taken together these suggest an OeBe classification. 
The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the pres¬ 
ence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally 
explained under the assumption that the accretor is a neutron star. 

Conclusions. VETS 399 appears to be the first high-mass X-ray binary identified within 30 Dor, sharing many observational charac¬ 
teristics with classical Be X-ray binaries. Comparison of the current properties of VETS 399 to binary-evolution models suggests a 
progenitor mass > 25Mo for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. 
VETS 399 is now the second member of the cohort of rapidly rotating ‘single’ O-type stars in 30 Dor to show evidence of binary 
interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar 
population. 

Key words, stars:evolution - stars:early type - stars:individual:VFTS 399 


1. Introduction 

Located within the Large Magellanic Cloud (LMC), the 30 
Doradus complex is the most luminous H ii region in the Local 
Group. It contains > 10^ OB stars, with star formation appar¬ 
ently commencing ~25Myr ago and continuing to the present 
da y (Wa lborn & Blades IT%7l Grebel & Chu l2^ Walborn et 
al. 120131) . Combined with a moderate line-of-sight extinction, 30 
Dor therefore provides an exceptional opportunity to investigate 
the complete life-cycle of massive stars. In order to exploit this 
potential, a multi-epoch spectroscopic observing campaign - the 
VLT-FLAMES Tarantula Survey (VETS) - was undertaken be- 

* Based on observations collected at the European Southern 
Observatory under program ID 182.D-0222 


tween 2008 October-2009 February; a full description of this 
investigation, including target list, observational strategy and 
data-reduction methods employed was presented by Evans et al. 

(ISOTTTl . 

Fulfilling a design goal of the programme, the multi-epoch 
observations enabled the identification of both single and binary 
OB star populations (Sana et al. 120131) . Investigation of the single 
O-type stars revealed that the distribution of their rotational ve¬ 
locities has a two-component structure (Ramirez-Agudelo et al. 
120131 ). comprising a low velocity peak (~ 80kms ') and a high- 
velocity tail (extending to ~ 600kms '). The origin of the high 
velocity cohort is currently uncertain, with the relative contribu¬ 
tion from natal and binary-driven channels still ill-constrained 
(de Mink et al. 120131 Ramirez-Agudelo et al. 120131) . With a pro- 
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Fig. 1. WFI image of the southern reaches of 30 Dor with 
VFTS 399 indicated, highlighting its location on the periphery 
of the complex. The held of view is ~ 135/7cxl25pc for a LMC 
distance modulus of 18.5mag (e.g. Gibson l2000t . 

jected equatorial rotational velocity, VgSin/ ~ 334 + 18kms“' 
(Ramirez-Agudelo et al. I2013I I. the 09 Illn star VFTS 399 is a 
member of this group (Walborn et al. 120141 1. while its location 
on the periphery of the 30 Dor complex (Fig. 1) is suggestive 
of a potential runaway nature (cf. Blaauw [l961l) . Unexpectedly, 
cross correlation of the VFTS sample with the X-ray survey of 
Townsley et al. (120141) revealed it to be both highly variable and, 
at its peak, one of the most X-ray luminous stars within 30 Dor. 
Motivated by these observational hndings we report the results 
of a multi-wavelength analysis of VFTS 399 in order to clarify 
its nature, presenting the data in Sect. 2, utilising these to infer its 
physical properties in Sect. 3, discuss the implications of these 
hndings in Sect. 4, before summarising our results in Sect. 5. 

2. The observational properties of VFTS 399 

2.1. Optical and IR properties 

2.1.1. Photometry 

Evans et al. (I201 It provided single-epoch optical and near-IR 
photometry of VFTS 399, which we present in Table 1. The 
latter are derived from the survey of the LMC by Kato et al. 
(12007b . which were obtained on 2003 November 16. Earlier 
2MASS magnitudes, obtained on 1998 March 19, are also avail¬ 
able (Skrutskie et al. l2006lF] and are suggestive of variability (see 
below). Einally Spitzer SAGE observations (Meixner et al. 12006b 
obtained between 2005 July 15-26 and October 26-November 2 
are also presented in Table 1. 

VETS 399 also lies within helds covered by the Optical 
Gravitational Lensing Experiment (OGLE-III; Udalski 12003b 
and VISTA-VMC (Cioni et al. 1201 lb surveys. Unfortunately, 
neither dataset coincides with the period over which spectro¬ 
scopic observations were made but both clearly indicate that 
VETS 399 is variable at the level of ~ 0.1 - 0.2mag over 
timescales of ~ 10 - 10^ days in the V, I and Xj-bands (Eig. 2). 

* J = 15.48, H = 15.61 and = 14.83, noting the Xj-band magni¬ 
tude is flagged as uncertain. 


Table 1. Optical and IR photometry of VETS 399. 


B 

V 

J 

H 

Xs 

[3.6] 

[4.5] 

15.91 

15.83 

15.40 

15.30 

15.22 

14.95 

14.95 


In the extensive /-band dataset the variability appears charac¬ 
terised by signihcant scatter superimposed on longer-term vari¬ 
ability. The Xs-band data are more sparsely sampled, but are 
suggestive of coherent variability over a ~ 10^ day timescale. 
Preliminary analyses of the more extensive OGLE-III /-band 
dataset revealed a number of potential periodicities - e.g. ~ 
156,370 and 720 days - although the later two appeared to result 
from the seasonality of the observations. Neveretheless, in or¬ 
der to assess the signihcance of all the periodicities we followed 
the methodology employed by Bird et al. (I2012|FI . No coherent 
statistically-signihcant periods were identihed for the complete 
dataset, while searches over subsets of these data also returned 
no convincing evidence for transient periodicity. 

2.1.2. Spectroscopy 

The VLT/ELAMES spectroscopy of VETS 399 is summarised 
in Table 2 and the data-reduction techniques employed are fully 
described in Evans et al. (1201 lb . noting that sky subtraction was 
accomplished utilising a median sky spectrum constructed from 
all the sky hbres of the frame that were uncontaminated by flux 
from bright stars. Walborn et al. (120141 their Eig. 10) presented 
the 3960-4750A portion of the spectrum from which they as¬ 
signed an 09 Illn classification; for brevity we do not replicate 
this here, but do plot the Ho', H/1 and He i 6678A emission line 
profiles in Eig. 3. 

All three lines demonstrate complex, morphologically simi¬ 
lar profiles with pronounced high-velocity emission, of compa¬ 
rable intensity in both blue and red line wings and a strong, nar¬ 
row component centred on the systemic velocity of VETS 399. 
Similar, single-peaked emission at low projected velocities is 
also present in the photospheric profiles of the higher Balmer 
series and select He i lines (Walborn et al. 12014b . Given the 
presence of spatially-variable nebular emission across 30 Dor 
we caution against the over interpretation of the profiles of the 
core components in these lines. Sky subtraction successfully re¬ 
moved the [N ii] 6547A nebular line, but left small single-peaked 
residuals in the [Nii] 6584A and [Sii] 6717-)-6730A lines. 
Consequently we suspect that nebular contamination might also 
be present at low projected velocities in the Ha, H/3 and Her 
6678A transitions. We discuss the physical interpretation of 
these profiles further in Sect. 3. 

Ramlrez-Agudelo et al. (12013b determined VgSin/ ~ 334 + 
18kms ' from the LR02 and LR03 spectra, while Sana et al. 
(12013b utilised the multi-epoch nature of the dataset to search for 
radial velocity (RV) shifts. They found that VETS 399 demon¬ 
strated epoch-to-epoch RV shifts with a peak-to-peak amplitude 
of ~ 50kms '. Given the magnitudes of the observational uncer¬ 
tainties (Table 2), this is significant at the ~ 2.7cr level, meaning 


^ A fast implementation of the Lomb-Scargle periodogram, with peri¬ 
ods searched for ranging from 2 days (the Nyquist frequency) to '-1000 
days (dictated by the length of the lightcurve). A rolling-mean was fur¬ 
ther calculated for the lightcurve and was used to de-trend the data, with 
extensive Monte-Carlo simulations employed to determine the signif¬ 
icance and associated errors of any periods detected. See Bird et al. 
l l2012b for further details. 


2 














































J. S. Clark et al.: The VLT-FLAMES Tarantula survey XX. The nature of the X-ray bright emission line star VFTS 399 



15.60 


15.55 


15.50 



15.85 


15.90 - 

2000 


5000 4000 

HJD-2450000 


5000 


Fig. 2. VISTA-VMC /Ts-band (upper panel) and OGLE 
/-(middle panel) and V-band (lower panel) lightcurves for 
VFTS 399. 


Table 2. Summary of spectroscopic observations. 


Mode - 1 - 
Setting 

/(-coverage 

(A) 

R 

Date 

RV 

(km s“') 

Medusa LR02 

3960-4564 

7000 

22/12/08 

266.8 ± 13.3 




24/12/08 

286.5 ±9.1 




27/01/09 

300.3 ± 17.3 




11 mm 

265.8 ± 16.0 




07/10/09 

276.5 ± 19.2 

Medusa LR03 

4499-5071 

8500 

20/12/08 

304.9 ± 10.6 




21/12/08 

317.1 ± 12.7 




22/12/08 

297.3 ± 13.4 

Medusa HR 15N 

6442-6817 

16000 

18/12/08 

- 




19/12/08 

- 


The He It 4200A, Hei 4387A and Hen 4541A transitions were em¬ 
ployed for RV determination with Medusa LR02 observations, and He i 
4713A and Hei 4922A with Medusa LR03. Note the LR02 and LR03 
observations made on 2008 Dec. 22 were separated by ~ 68mins. 


that VFTS 399 fell below the conservative threshold set by these 
authors for classihcation as a binary. We return to this issue in 
Sect. 3.1, noting that these observations also yield a mean RV 
value of 289 + 19kms“'. 



Fig. 3. Montage of the Ha, fi and He i 6678A line prohles, with 
the latter two plotted with an arbitrary vertical exaggeration and 
offset in order to ease comparison to Ha. Note the central com¬ 
ponent of Ha at projected velocity ~ 0km s ' has a maximum 
peak intensity of ~ 6x continuum; we suspect nebular contam¬ 
ination of this feature and the comparable peaks in H/? and He i 
6678A. 


2.2. X-ray properties 

Fortuitously, VFTS 399 has been the subject of a number of 
X-ray observations spanning ~ ISyears. Both XMM-Newton 
(XMMU J053833.9-691157; Shtykovskiy & Gilfanov l2005t and 
Chan dra (C XOU J053833.4-691158 = [TBF2006] 27; Townsley 
et al. I2006t provide X-ray point source detections that may be 
associated with VFTS 399. Utilising data from 2001 November 
19, the former study quotes L 2 -iokeV ~ LI x lO^^erg s“', as¬ 
suming a power law spectrum with photon index, T = 1.7 and 
Ah = 6 X lO^^cm^^. The latter study utilised a ~ 21ks observa¬ 
tion from 1999 September and returned Lo.5-8keV ~ 10^"'erg s“' 
following a power law ht with T =1.7 and Ah = 1.3 x 10^^cm“^. 

Subsequently, further Chandra observations were made on 
2006 January 21, 22 and 30 (Townsley et al. l20T4ll . VFTS 399 
is also located within the target held of the current ~ 2Ms T- 
ReX Chandra programme (PI; L. Townsley) which runs through 
2015; consequently we appraised the multiple additional obser¬ 
vations available to us at the time of writing (Table 3), deferring 
an analysis of the full data set for a future paper. 

VFTS 399 was detected in each observation and was found 
to be well separated from other X-ray sources. The resultant 
lightcurve is presented in Fig. 4 and it is immediately apparent 
that the source is highly variable, even within a single obser¬ 
vational segment. In order to screen for variability in the gross 
morphological properties of the spectra we determined the me¬ 
dian event energy for individual observations (Fig. 4). This met¬ 
ric shows no evidence for spectral variability across this ob¬ 
servational sequence. In order to consolidate this hnding we 
analysed the hrst eight epochs of data which, when taken to¬ 
gether, appear representative of the whole dataset inasmuch as 
they fully sample the full dynamical range found for VFTS 399 
(Lo.5-8keV ~ 10^''-10^^erg s“'). Furthermore, excludingthe 1999 
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Table 3. Summary of Chandra X-ray observations of VFTS 399 
during 2006 and 2014. 


Obs. ID 

Start Time 

Exposure (s) 

5906 

2006-01-21T19:04 

12317 

7263 

2006-01-22T16:51 

42528 

7264 

2006-01-30T15:06 

37593 

16192 

2014-05-03T04:10 

93760 

16193 

2014-05-08T10:15 

75993 

16612 

2014-05-11T02:15 

22672 

16194 

2014-05-12T20:00 

31335 

16615 

2014-05-15T08:24 

45170 

16195 

2014-05-24T14:09 

44405 

16196 

2014-05-30T00:05 

67109 

16617 

2014-05-31T01:27 

58860 

16616 

2014-06-03T22:26 

34529 

16197 

2014-06-06T12:32 

67802 

16198 

2014-06-11T20:20 

39465 

16621 

2014-06-14T14:46 

44399 

16200 

2014-06-26T20:01 

27360 

16201 

2014-07-21T22:13 

58393 

16640 

2014-07-24Tll:21 

61679 

16202 

2014-08-19T15:30 

65128 

17312 

2014-08-22T06:21 

44895 


data due to calibration issues, we also summed all remaining 
observations to produce a single, high S/N spectrum which we 
treated in the same manner. 

Specifically, we first subtracted a background spectrum (cor¬ 
responding to ~ 3% of the extracted 0.5-8.0keV X-ray counts) 
before several different models within xspec (Arnaud 119961 1 
were employed in an attempt to fit these data. A simple absorbed 
power-law (tbvarabs*pow) with Z = 0.4 x Zq was found to 
provide the best fit (Fig. 5): parameters returned for this model 
upon application to the summed spectrum were F = l O^o p 
Nh = 2.7;^q3 X lO^^cm^^ (corTesponding to Ay ~ 16.9^j gmag, 
suggestive of signihcant obscuration intrinsic to VFTS 399) and 
an absorption corrected Lq.s- Skev ~ 5.0 X lO^'^erg s '. 

The results obtained from analysis of the first eight indi¬ 
vidual epochs are presented in Table 4; this yielded a range 
of parameters spanning 0.8 < F < 1.1, 2.4 x 10^^ < Nu < 
3.1 X lO^^cm^^ and 1.0 x 10^"^ < To.s-skev ^ 12.6 x 10^"^erg s“*. 
For epochs with very few counts Nh and F were held at the value 
determined from the summed spectrum, with only the normalisa¬ 
tion allowed to vary. We caution that fits are provided for epochs 
as a whole and hence do not account for the variable behaviour 
of the source over the course of a given observation (cf. Fig. 4). 
Moreover, while some datasets are htted by this simple model, 
others are not; this likely indicates dehciencies in the model, but 
the low count rate precludes more sophisticated treatments at 
this stage. Nevertheless, we conclude that there is no evidence 
for evolution of the gross morphological properties of the X-ray 
spectrum of VFTS 399 over the course of the observations. 

2.2.1. X-ray timing analysis 

Lomb-Scargle (LS; Lomb|T97^ Scargle[T982l Press & Rybicki 
I1989b analysis was performed on each of the single observation 
light curves as well as on a composite light curve, made up of 
all the data sets. LS analysis is ideal for detecting weak, periodic 
modulation in unevenly sampled data sets (such as our compos¬ 
ite data set). A ’blind’ search for periods between 7s (approx¬ 
imately twice the frame interval of Chandra) and IxlO^s was 



Fig. 5. Top panekAbsorbed power-law fit to the summed spec¬ 
trum constructed from the 20 epochs of Chandra observations 
shown in Table 3 (5397 net counts; 975ks of combined expo¬ 
sure). Bottom panel; residuals for this fit. 



Fig. 6. LS periodogram for the combined dataset, the dashed line 
representing the Act confidence interval, which the peak corre¬ 
sponding to the ~ 2567 s period clearly exceeds. 


performed with a frequency step of 1 x 10“^. Foreshadowing the 
discussion in Sects. 3.2 and 4.1, this exceeds the observed range 
of pulsational periods exhibited by neutron stars in X-ray bina¬ 
ries with OB star primaries. 

At the other extreme, long-term modulation of the lightcurve 
might be expected from the rotation of magnetic stars, or due to 
orbital modulation in colliding-wind systems or accreting high- 
mass X-ray binaries; in all cases the resultant period may range 
from several days to several years. Unfortunately, the durations 
of the individual observing blocks are < Iday (with only a sin¬ 
gle exception; Table 3) meaning that none would be expected 
to sample a full cycle of a system with a period of a few days. 
Moreover, the overall time span of the observations is such that 
we would also not expect to fully sample a single cycle of any 
long-period (> lO^day) modulation. Thus, we conclude that the 
current dataset is ill-suited to the robust detection of any such 
putative period; a longer (continuous) time base of observations 
would be required. 

Furthermore, orbital modulation of the X-ray light curves 
of such binaries occurs in the subset of eccentric-orbit systems. 
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Time (ks) 


Fig. 4. Concatenated Chandra observations of VFTS 399 (975 ks combined exposure). Binned flux and associated Icr errors are 
shown in black, while binned median energies are shown in red. A sliding-window lightcurve inferred from the individual events 
(i.e. no binning) is shown by the blue lines, with a 68% pointwise confidence band limit represented by the two lines. Running left 
to right, the start dates of the 20 epochs of observations are given in Table 3. 


which may demonstrate high-amplitude flaring during perisatron 
passage due to strengthened wind/wind interactions in colliding 
wind systems (e.g. Corcoran I2005t or enhanced accretion rates 
in X-ray binaries (e.g. Galache et al. I2008t . Such repetitive, pe¬ 
riodic flaring is not immediately apparent in the lightcurve of 
VFTS 399, although this does not preclude its future occurrence 
since, in the case of X-ray binaries, such behaviour is observed 
to be transient. 

With the above constraints in mind. Fig. 6 shows the LS pe- 
riodogram of the whole data set. There appears to be evidence 
for a period at ~ 2567 s at > 4cr level of significance. Fig. 7 
shows the pulse profile, which appears asymmetric and single 
peaked. A peak at this period is present in the periodograms of 
a number of individual observation blocks, albeit inevitably at 
a lower level of significance. The strength of the detection of a 
period in a given observation depends on the S/N of the obser¬ 
vation, the pulsed fraction and the period length: it is easier to 
And shorter rather than longer periods in any given observation 
since more ’cycles’ are present. As such it is not surprising that 
this period is not robustly detected in some observations; for in¬ 
stance our shortest observation, 5906, is only 12.3 ks in duration. 


which allows for <5 pulsations. Nevertheless, this latter result 
gives us confidence that the period is not spuriously introduced 
by our analysis of the composite lightcurve during the process 
of ‘stitching together’ its constituent parts; this assumption was 
quantitatively tested as described belo'wQ Finally, a strong pe¬ 
riod at 707 s detected in the periodogram of observation 16193 
is the Chandra dither period, resulting from a bad CCD column 
falling within the aperture. 

In order to test both the robustness and statistical significance 
of this period a number of simulations of the data were run. LS 
analysis was first performed on the barycentric-corrected X-ray 
event timestamps with a variety of sizes to make sure that the 
2567 s period was not an artefact of the binning adopted - in 
each case the period was still detected. The data were then boot¬ 
strapped to determine how robust the period was. This was ac- 

^ Such a methodology is routinely applied in e.g. the case of optical 
observations of X-ray binaries, where periodic gaps between distinct 
segments of the lightcurve are introduced by seasonal visibility con¬ 
straints (e.g. Bird et al. 12012b . Moreover, Reig et al. 12009 b encountered 
exactly this issue in their search for X-ray pulsations in the high-mass 
binary 4U2206+54, due to visibility gaps in their RXTE lightcurve. 
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Table 4. Summary of indicative fits to the Chandra X-ray observations of 2006 and 2014. 


ObsID 

Date 

Exposure 

Net 

Ah 

F 

Normalisation (x 10 = 

T^/DoF 

logTfl.S-Skev 



(s) 

Counts 

(xlO^^ cm-2) 


photons cm“^ s“' keV“‘) 


(ergs^') 

Summed 

- 

975400 

5397 

9 7+0.4 

03 

1.0-i 

1.5— ■“ 

125.5/115 

34.7 

5906 

21/01/06 

12320 

175 

1.5 

0.8-;= 

9 q+3.0 

4.8/8 

35.1 

7263 

22/01/06 

42528 

407 

1.0 

1-oil 

2 4+1.4 

0.8 

24.1/23 

34.9 

7264 

30/01/06 

37593 

40 

2.7 

1.0 

0 2“^^ * 

0,1 

6.7/9 

34.0 

16192 

03/05/14 

93760 

419 

1.3 

03 

1 

0 + 1.9 

29.6/24 

34.7 

16193 

08/05/14 

75996 

657 

-1.0 

1 

^•^-0.3 

38.6/38 

35.0 

16612 

11/05/14 

22678 

54 

2.7 

1.0 

0-6 02 

1.8/ 

9.8/7 

34.4 

16194 

12/05/14 

31335 

238 

2.7 

OKI 

30.6/13 

34.8 

16615 

15/05/14 

45170 

237 

9 4+1-8 

-0.5 

6.4/12 

34.7 


Summary of the X-ray model parameters for the summed spectrum (top row) and the individual spectra derived from the first eight representative, 
constituent spectra. Parameters given in italics have been fixed to the properties derived from the summed spectrum due to low count rate, note 
also the comparatively poor fit to the spectrum of 2014 May 12. In all cases errors quoted are 90% confidence limifs. The X-ray luminosities are 
derived after correction for extinction and for a distance to the LMC of SOkpc. 



Phase 

Fig. 7. X-ray data folded on the ~ 2567 s period. The pulsed 
fraction, calculated by integrating over the profile and computing 
the ratio between the pulsed section of the profile and the total 
flux, is 0.17 + 0.04. 


complished by creating ‘fake’ data sets by using a random num¬ 
ber generator to select ~ 63% of the data points and performing 
a LS analysis on the resultant lightcurve. If the detected period 
was being driven by a small number of data points it will not 
be consistently detected during the bootstrapping. The period at 
~ 2567 s was detected in 96.6% of the trials, suggesting that the 
period is indeed robust. 

We also ran a second bootstrapping analysis, this time select¬ 
ing random observations (as opposed to data points) to build the 
simulated data sets. Anywhere from 30 to 80% of the total data 
were used (since there are fewer ways to randomly select 20 light 
curves in contrast to the full set of data points). In this instance 
the period at ~ 2567 s was recovered in 85% of the trials. The 
discrepancy between the two analyses almost certainly arises in 
the low data fraction regime in the second bootstrapping, since 
recovering the period with only 30% of the original data will be 
unlikely. 

The significance of the detected period was assessed by esti¬ 
mating the probability of the detection under the null hypothesis 
that the lightcurve is uniform. Over 15000 such light curves were 
generated by randomizing the X-ray event timestamps. This al¬ 


lows us to create simulated light curves with the same statistical 
properties as the original data set. LS analysis was performed on 
each of these light curves and the maximum LS power achieved 
recorded. This is shown as the dashed line in Fig. 6. The re¬ 
sults of the Monte-Carlo analysis suggest that our ~ 2567 s pe¬ 
riod, with a LS power of 42.8 has a less than 1 in 15787 (i.e. 
<0.006334%) probability of being a chance occurrence due to 
white noise. Although a second peak (at 8718 s) in the peri- 
odogram exceeds this significance threshold, it is not robust to 
the above bootstrapping analyses; specifically, the peak disap¬ 
pears into the noise if the three observing blocks from 2006 are 
excluded, suggesting that it is most likely an artefact. 

We therefore conclude that the Monte-Carlo and bootstrap¬ 
ping analyses both suggest there is a periodic modulation in our 
data set although, trivially, it cannot reveal the origin of this peri¬ 
odicity. The ~ 2567 s period does not correspond to the Chandra 
frame interval (3.241040 s), dither period (707 s) or the orbital 
period (198 ks), nor does it correspond to an alias (the beat fre¬ 
quency) of any of these three potential periods with each other 
(which corresponds to 3.241093, 3.255966 s and 707.5335 s; for 
a comprehensive discussion on aliasing and periods see Bird et 
al. l2012l l. For the remainder of the paper we therefore proceed 
under the assumption that the period is astrophysical in origin. 

Typically, in determining the error on a period one assumes 
that the period derivative is zero throughout the observation and 
hence that the dominant source of errors is the Poisson error on 
the counts themselvefl However, given the time span of our ob¬ 
servations we may not a priori make this assumption. This in 
turn prevents a robust quantitative estimate of the error on the 
period since we may not distinguish between random error and 
a systematic evolution of the period during the course of the ob¬ 
servation. 

However, this turns out to be of secondary concern as it is 
possible that the ~ 2567s period we detect could be an alias 
of the true astrophysical period and the Chandra frame inter¬ 
val (3.241040 s; note that aliasing is a common problem in the 
analysis of e.g. optical light curves of X-ray binaries (Bird et al. 
120121) ). This would require the true period of VFTS 399 to be 

In such circumstances one may vary the lightcurve within the errors 
on the individual data points in order to produce a large number (i.e. ~ 
10^) of synthetic lightcurves that are statistically identical to the original 
and subject each to a LS search. A histogram of the resultant distribution 
of periods is typically well approximated by a Gaussian, from which 
the width (i.e. standard deviation) yields a quantitative estimate of the 
period error. 
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~ 3.24515 s. Unfortunately following the above discussion such 
a period would not be directly detectable in our curi'ent dataset; 
observations by a different observatory such as XMM would be 
necessary to distinguish between these possibilities. It is, how¬ 
ever, noteworthy that our simulations indicate the ~ 2567s pe¬ 
riod is robust to evolution of the period, such as those commonly 
seen in X-ray binaries due to accretion torques and doppler mo¬ 
tion, but not so if it is due to an alias. Given the nature of the 
system (Sect. 3.2 and 4.1) this might be taken as evidence to 
disfavour the alias hypothesis. 

Nevertheless, it is important to emphasise that either inter¬ 
pretation still requires a real periodicity in the X-ray lightcurve 
to be present. However, we are left in the uncomfortable posi¬ 
tion of having identihed a statistically-signihcant period of ap¬ 
parently astrophysical origin in our dataset but not being able 
to conclusively determine whether the ~ 2567 s periodicity is 
an alias or not. As a consequence for the remainder of the pa¬ 
per (e.g. Sects. 3.2 and 4) we are forced to consider the physical 
implications of both the ‘short’ and ‘long’ periods in parallel. 
Clearly, independent observations of VFTS 399 with observato¬ 
ries other than Chandra would be highly desirable, since one 
requires a much shorter frame time in order to test the alias hy¬ 
pothesis. 


3. The nature of VFTS339 

VFTS 399 is displaced from the locus occupied by the majority 
of OB stars in both near-IR {{J - H)f(H - Xj); Dougherty et 
al. lT994ll and mid-IR ([3.6]/([3.6]-[4.5]); Bonanos et al. l2^ 
colour/colour and colour/magnitude diagrams. It is in a region 
populated by, amongst others, early-type stars exhibiting the 
classical Be phenomenon. Furthermore, photometric variability, 
such as that demonstrated by VFTS 399, is a well-known prop¬ 
erty of classical OeBe stars (e.g. Mennickent et al. 1200211 . 

Emission in the Balmer series is another key classification 
criterion for classical OeBe stars. Under such a scenario, the 
emission at large projected velocities in the line wings of Ha, 
Hy6 and Hei 6678A (Fig. 3) would result from the presence 
of a compact, quasi-Keplerian decretion disc; during the early 
stages of disc-reformation, time resolved Ha spectra of o And 
revealed emission components with a comparable peak-to-peak 
separation (Clark et al. I2003t . Similar behaviour is also ob¬ 
served in both the Ha and Hei 6678A lines of the BO Ve pri¬ 
mary of the high mass X-ray binary X Per. Comparison of X 
Per to VFTS 399 appears particularly apposite, given that both 
stars demonstrate additional emission features interior to those 
present in the line wings. In X Per this is attributed to the 
presence of a gaseous circumstellar ring exterior to, and rotat¬ 
ing more slowly than, the inner disc that is responsible for the 
high-velocity emission component of the line prohle (Tarasov & 
Roche [T9951 Clark et al. lTODTlR 

The OeBe phenomenon extends to stars with spectral types 
as early as that of VFTS399 (e.g. Neguerueula et al. 120041 1. while 
such stars are also known to be rapid rotators (e.g. Townsend et 
al. 120041 1. In light of these observational properties, VFTS 399 


^ In this regard we are mindful of the likelihood of nebular contam¬ 
ination in VFTS 399. However even if all emission interior to the high 
velocity wings were nebular in origin - which appears unlikely upon 
comparison of the relative widths of the nebular lines to this component 
- the global consideration that VFTS 399 is a classical OeBe star would 
not be affected. However, potential unresolved nebular contamination 
does preclude the measurement of the line equivalent widths. 


appears to be a classical OeBe star and hence we slightly revise 
its formal classihcation to 09 Illne. 


3.1. Stellar parameters 

The broad-band photometry presented in Table 1 was anal¬ 
ysed using the Bayesian code CHORIZOS in order to deter¬ 
mine the interstellar reddening (Malz Apellaniz 120041 Mafz 
Apellaniz et al. 120141 1. from which Walborn et al. (120141) infetTed 
My ~ -4.41, Teff ~ 32.8kK and \og{L\,o\/L@) ~ 4.9 + 0.1 for 
VFTS 399. As part of a wider programme (Ramfrez-Agudelo, 
in prep.), the LR02 and LR03 spectroscopic data were subse¬ 
quently analysed with the non-local thermodynamic equilibrium 
code FASTWIND (Puls et al. [20051) . yielding ~ 31.2+^ °kK, 
logg ~ 3.6and log(Lboi/To) ~ 4.8 + 0.1. 

We have also attempted to estimate the nitrogen abundance 
in VFTS 399. Due to its large projected rotational velocity, the 
N II spectrum could not be analysed but the N iii multiplet near 
4515A was observed as a blend in the LR02 and LR03 spectra. 
In both spectra, this region lies close to either the upper or lower 
wavelength limits making normalisation of the spectra difficult. 
We have therefore chosen to treat the spectra separately, with 
a simple subjective normalisation. These spectra are shown in 
Fig. 8, together with predictions from model atmosphere calcu¬ 
lations. The latter were taken from the O-type star grid of Lanz 
& Hubeny (120031 ) and have been convolved with an appropriate 
broadening function. The weaker theoretical N iii spectrum has 
been interpolated to the FASTWIND atmospheric parameters 
(and a microturbulence of lOkms ') and is for a nitrogen abun¬ 
dance of 7.62 dex and a EMC metallicity. The stronger spectrum 
is for a nitrogen abundance of 7.92 dex but a galactic metallicity; 
this leads to more line blanketing, a hotter temperature structure 
and, for example, enhanced He ii absorption. Inspection of grid 
models with different effective temperatures implies that this will 
lead to an overestimation of the strength of the N iii spectra com¬ 
pared with an EMC metallicity. 

Both the LR02 and LR03 spectra imply a nitrogen abundance 
of approximately 7.8 dex. However we emphasise that this value 
and the atmospheric parameters from the FASTWIND analysis 
must be treated with caution due to the limitations of classical 
non-LTE techniques. Both the Ha and H/3 profiles are subject 
to inhlling, with some evidence for similar behaviour in the Hy 
prohle, although a dehnitive conclusion is complicated by pos¬ 
sible nebular contamination. Moreover we have assumed that 
an isolated stellar model atmosphere is appropriate; the possi¬ 
ble systematic errors introduced by the complexity of the system 
e.g the presence of an accreting secondary (see Sect. 3.2) are 
difficult to quantify but could be signihcant. Nevertheless, we 
are conhdent that our conclusion that VETS 399 demonstrates a 
pronounced nitrogen enhancement is robust. 

Likewise, contamination of the broad-band photometry by 
continuum emission from the circumstellar disc inferred for 
VETS 399 has the potential to seriously compromise reddening, 
and hence luminosity, estimates. This is evident upon consid¬ 
eration of the long term photometric behaviour of the BeX-ray 
binary A0535-I-26. Over a ~ 16yr period Clark et al. (119991 ) re¬ 
port Ay ~ 0.4mag and AE{B -V)~ 0.13mag as a result of disc 
variability (in the sense the system became redder as it bright¬ 
ened). As a consequence if we fail to account for the presence of 
the circumstellar disc, we risk overestimating the bolometric lu¬ 
minosity of VETS399 because (i) we will not subtract the excess 
disc continuum emission and (ii) because we will overestimate 
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the interstellar reddening component and hence over correct for 
this. _ 

Riquelme et al. (120121) investigated this problem for Be X- 
ray binaries in the Galaxy and Magellanic Clouds, hnding it to 
be signihcant for the latter objects, where contributions from 
circumstellar and interstellar reddening are comparable. In the 
absence of contemporaneous spectroscopic and photometric ob¬ 
servations of VFTS 399 in a discless state we may not determine 
the interstellar reddening. Moreover, nebular contamination of 
the Ha line prevents an equivalent width measurement and from 
this an estimate of the circumstellar reddening via the relation¬ 
ship dehned by Riquelme et al. (120121 ). As a consequence we 
are forced to adopt a ‘worst-case’ approach, correcting for the 
maximum contribution from circumstellar reddening reported by 
Riquelme et al. (120121 ) of - V) ~ 0.17mag. Utilising this 
value, the photometry from Table 1 and the intrinsic colours of 
an 09 III star (Martins & Fez I2006t we may infer an interstel¬ 
lar reddening. Then employing the relation from Riquelme et al. 
(12012b between total optical extinction and interstellar redden¬ 
ing - = 3.1E‘%B -V)- 0.3mag - for Be X-ray binaries and 

an appropriate bolometric correction (from Martins et al. 12005b 
we determine log(Lboi/To) ~ 4.4 for VFTS 399. 

This value appears surprisingly low in comparison to either 
log(Lboi/To) ~ 5.17or4.72 expected for, respectively, 09 giants 
and main sequence stars (Martins et al. 12005b . It is difficult to 
envisage any physical process, such as binary interaction, that 
would leave a star this under luminous for its spectral classih- 
cation, suggesting that we have over corrected for circumstel¬ 
lar contamination. Therefore we favour the intrinsic luminosity 
of an 09 V star as an appropriate lower bound for VFTS 39S0 
and the upper bound from Walborn et al. (12014b . leading to 
log(Lboi/To)~4.7-4.9. 

We employed the Bayesian code BONNSAfl (Schneider et 
al. I2014bb to determine the current mass (Mcurr) and age of 
VFTS 399 (Table 5). Firstly, we neglected the luminosity esti¬ 
mate and simultaneously matched the observed effective tem¬ 
perature, Teff = 31.2;|;p°kK, surface gravity, logg = 3.60j;Q[g and 
projected surface rotational velocity, Ve sin / = 334 + 33kms“' to 
the rotating, single star LMC models of Brott et al. (1201 lb and 
Kohler et al. (12014b . We used a Salpeter initial mass function 
(Salpeter ll955b and the observed distribution of rotational ve¬ 
locities of the 30 Doradus O-type stars (Ramirez-Agudelo et al. 
12013b as appropriate priors, together with a uniform distribution 
of stellar ages and a random orientation of rotational axes. The 
stellar models predicted a luminosity of log(Lboi/To) = 5.0;|;gjg, 
in agreement with the luminosity estimate of Walborn et al. 
(12014b and the expectations from the 09 classihcation. A lumi¬ 
nosity of log(Lboi/To) ~ 4.4 was excluded at more than 3cr. 

Secondly, we took the luminosity estimates for VFTS 399 
into account. Following from above we assumed log(Lboi/To) ~ 
4.7 - 4.9 but, for completeness, also undertake calculations as¬ 
suming both log(Lboi/To) ~ 4.4 - appropriate for the max¬ 
imum expected contamination by a circumstellar disc - and 
log(Lboi/To) ~ 5.2 - the expected luminosity for an 09 III star 

® Foreshadowing the discussion in Sect. 4.1, Negueruela & Reig 
<20011 reported that the emission-line mass-donor BD+53“2790 (clas¬ 
sified as 09.5 V) in the high-mass X-ray binary 4U 2206-1-54 demon¬ 
strates variable emission, infilling the He ii 4686A photospheric line that 
forms part of the luminosity classification criteria for late O-type stars 
(Walborn et al. 12014b . If such an effect is also present here then it could 
erroneously lead to the assignment of a giant rather than dwarf lumi¬ 
nosity class for VFTS 399. 

’ The BONNSAI web-service is available at 
http://www.astro.uni-bonn.de/stars/bonnsai 
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Fig. 8. LR02 (upper) and LR03 (lower) spectra for VFTS 399 
encompassing the Niii multiplet near 4515A. Also shown are 
theoretical spectra with nitrogen abundances of 7.62 dex (red, 
dash) and 7.92 (blue, dot-dash), convolved with an appropriate 
broadening function. The absorption near 4541A is due to He ii. 

(Martins et al. 12005b . The results are presented in Table 5. For 
log(Lboi/To) ~ 4.4, the stellar models could not reproduce the 
observables simultaneously, likely due to the low luminosity. For 
our favoured range of bolometric luminosities we determined 
Afcurr ~ 18.0 - 20.0/Wo. Finally, we inferred a current age of 
~ 6Myr for VFTS 399, assuming single-star evolution. However, 
pre-empting the discussion in Sect. 3.2, we consider it likely that 
VFTS 399 has accreted mass from a binary companion in the 
past, rejuvenating the star (e.g. Schneider et al. I2014ab : there¬ 
fore we consider this value as a lower limit for its true age. 

Under this assumption it is of interest to determine whether 
the observed nitrogen enhancement in VFTS 399 is the result 
of binary mass transfer. Langer (12012b showed that for Galactic 
stars such interaction may result in an enrichment of ~ 0.5 - 
O.Sdex. Since the initial C/N ratio is ~ 0.2 - 0.3dex larger in 
the Galaxy in comparison to the LMC we would expect the 
nitrogen enrichment in the LMC to be larger by this amount; 
for a baseline nitrogen abundance of ~ 6.9dex in the LMC we 
would predict a post-interaction abundance of ~ 7.6 - S.Odex. 
Unfortunately, while this is consistent with our findings for 
VFTS 399, it is also consistent with the predictions for rotational 
mixing in a single, rapidly rotating star, with BONNSAI predict¬ 
ing ~ 7.76 ^qj 2 - 7.83j/Qj5 dex for our favoured range of lumi¬ 
nosities. Therefore we are current unable to utilise the chemical 
abundance of VFTS 399 in isolation to distinguish between sin¬ 
gle and binary star evolution. 
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Table 5. Bonnsai predictions for the cutTent mass and lower limit 
for the age of VFTS 399 


log(Lbo,/Lo) 

Menrr 

Age 


(Wo) 

(Myr) 

4.7 ±0.1 

4.9 ±0.1 

5.2 ±0.1 

21.4+j-^ 

18.0?i| 

20.0+|« 

24.0?“ 

> 5.8?;;:? 

> 6 l+“ 

>6:if 


All four models utilise stellar temperature, surface gravity and pro¬ 
jected equatorial rotational velocity as input parameters, while rows 2-4 
present the results of including the stellar luminosity estimates (column 
1) as an additional constraint. The quoted errors are Icr uncertainties. 


3.2. Understanding the X-ray emission 

Single OB stars are known to be (thermal) X-ray emitters, likely 
due to shocks propagating through their winds. However, this 
emission scales as log(Lx/f'boi) ~ 10“^ (e.g. Berghoefer et al. 
119971 Sana et al. l2006t . while VFTS 399 is orders of magnitude 
brighter, with log(Lx/7'boi) ~ 0.3 - 4 x lO^'* * in the observations 
reported above, while the emission is of a non-thermal nature. 
These properties also exclude membership of the class of X-ray 
bright and variable y Cas-like OeBe stars, for which a high tem¬ 
perature thermal spectrum and log(Lx/7'boi) ~ 10“^ are expected 
(Smith et al. 120121 Rauw et al. 120131) . Likewise, the peak X- 
ray luminosity of VFTS 399 is two orders of magnitude greater 
than that found for Of?p stars and the subset of strongly mag¬ 
netic O-type stars which, in turn, are X-ray over-luminous with 
respect to the wider population of massive stars (cf. Naze et al. 
120081 Clark et al. |20()9| and refs, therein). Additionally, the X-ray 
emission in such stars is once again thermal in origin, in contrast 
to VFTS 399. 

At its peak luminosity, VFTS 399 is an order of magnitude 
brighter than any Oh-O colliding wind binary, although seven 
binaries with luminous blue variable or Wolf-Rayet primaries 
demonstrate logLx ^ lO^^erg s“^ (Gagne et al. l2012lFI . However 
in each of these systems the primaries support signihcantly more 
powerful winds than are expected for VFTS 399, suggesting that 
its X-ray flux is unlikely to arise via wind collisions (indeed no 
signature of any secondary is visible in our optical spectroscopic 
observations). Moreover the X-ray emission from such systems 
is thermal in origin, with kT ~ 0.9 - 4.4keV for those massive 
binaries with comparable X-ray fluxes. 

We are then left with the possibility that VFTS 399 contains 
an accreting relativistic object and, given the nature of the pri¬ 
mary, an obvious identification is with a Be X-ray binary. These 
systems are overwhelmingly comprised of OeBe stars (~09-B2 
V-III; Negueruela & Coe 120021 McBride et al. 120081) orbited by 
neutron stars and have periods 20 days < Porb ^ 400 days and 
eccentricities 0.1 < e < 0.9 (Cheng et al. 120141) . The interplay 
between orbital phase and the intrinsically variable circumstel- 
lar disc leads to complex, variable X-ray emission (e.g. Okazaki 
& Negueruela I200in Haberl et al. (120081 ) determined the X- 


WR48a, Mk34, rj Car, R140a, WR25, R136c and CXO J1745-28. 

* Most Be X-ray binaries are transient systems, with outbursts typi¬ 
cally > lOx their quiescent X-ray luminosities. Following the nomen¬ 
clature of Stella et al. ) ll986b . Type I outbursts occur periodically during 
periastron passage and typically yield fluxes of Lx ~ 10^’® - 10^^ erg 
s *. Type II outbursts are much rarer, aperiodic, have a longer duration 
and reach Lx > lO'^^erg s“*. In contrast, a subset of systems with wide 
(Torb ^ 100 days), low-eccentricity orbits are persistent sources with 
moderate luminosity (Lx ^ lO'^'* - 10^^^erg s“') and comparatively low- 


ray properties of a representative sample of SMC systems from 
single-epoch observations, finding the majority of spectra to be 
best fit by power laws (F ~ 0.93j;Q22) with a wide range of ab¬ 
sorbing columns (~ 10^°-10^^ cm“^) and luminosities spanning 
~ 1.5xl0^^-5.5xl03®ergs-'. The X-ray properties ofVFTS 399 
therefore appear entirely consistent with its classification as a Be 
X-ray binary. 

In such a scenario a neutron star accretor would be strongly 
favoured on both theoretical and observational grounds. Only 
one Ben-black hole X-ray binary has been identified from a 
Galactic population of ~80 such systems (Casares et al. 12014b . 
while the same study shows it to have a quiescent X-ray lumi¬ 
nosity some three orders of magnitude lower than VFTS 399. 

Mindful of the caveats in Sect. 2.2.1, the detection of pe¬ 
riodic modulation of the X-ray lightcurve provides support for 
such a conclusion, being most naturally explained as the pulsa- 
tional period of a rotating, accreting neutron star. Importantly, 
both ‘short’ (~ 3.245s) and ‘long’ (~ 2567s) potential period¬ 
icities lie within the range of pulsational periods exhibited by 
neutron stars within Be X-ray binaries (e.g. Knigge et al. l201 II 
Reig et al. 12009b and we expand upon this further in Sect. 4.1. 

4. VFTS 399 in context 

An immediate question to ask is whether we may constrain the 
orbital parameters of the system from the extant data? Adopting 
4Lcurr ~ 18.0 - 20.0Mq for the system primary (Sect. 3.1), the 
observed RV variability (Table 2 and Sect. 2.1.2) implies an un¬ 
expectedly short orbital period of Porb ~ 4 days for a canon¬ 
ical 1.4Mo neutron star companion. For Porb ~ 20days, con¬ 
sistent with the commonly found lower bound for Be X-ray bi¬ 
naries (Cheng et al. 12014b . we would infer a companion mass 
of ~ 2.5Mq, with any longer period strongly favouring a black 
hole, apparently excluded by the detection of X-ray pulsationJ^ 
Preliminary Monte-Carlo analysis following the methodology of 
Sana et al. (12013b appears to favour a lower-mass companion 
with a short period. We caution that this finding should be re¬ 
garded as provisional since it is primarily driven by the discrep¬ 
ancy in RV between two epochs separated by 0.8days (Table 2; 
2008 December 21, 22). Clearly further observations in order to 
build an orbital RV curve for the system are required to address 
this issue. 

4.1. The putative neutron star accretor 

Is the range of orbital periods suggested by the preceding anal¬ 
ysis consistent with one or both of the potential pulsational pe¬ 
riods? Following the relationship between the pulsational (Pspin) 
and orbital (Porb) periods in Be X-ray binaries delineated in the 
Corbet diagram (Corbet [l986l Knigge et al. l201 lb we would in¬ 
fer Port ^ 20 d for Pspin ~ 3.245 s and Porb ~ 500 days for 
Pspin ~ 2567 s. While the first combination appears internally 
consistent, the long orbital period implied by the second would 
appear to be in tension with our RV data. 

However, with Porb ^ 19.25 d and Pspin ~ 5560s (Reig et 
al. 12009112012b . the 09.5 Ve-i-neutron star binary 4U 2206-1-54 


level variability (< lOx quiescent flux), of which X Per is an exemplar 
(Haberl et al. fT^ Reig & Roche llOPTb . 

For comparison Antoniadis et al. (12013b determine ~ 2.01 ±0.04Mo 
for the pulsar PSR J0348-I-0432, while Clark et al. J2002b report ~ 
2.44 ± 0.21 Mq for the non-pulsating relativistic accretor in 4U1700-37 
and Casares et al. ( 12014b find ~ 3.8 - 6.9Mq for the black hole orbiting 
the Be star MWC656. 
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demonstrates that such a combination of long pulsation but short 
orbital period is viable, although this would require VFTS 399 to 
host the third most slowly rotating neutron star of any high mass 
X-ray binary, behind 4U 2206-H54 and 2S 01 14-H65 0 (Porb ~ 
11.59day and Pspin ~ 2.7hr; Crampton et al. 119851 Corbet et 
al. 119991 Farrell et al. 120081 1. In this regard we note that both 
the X-ray luminosity and aperiodic variability of VFTS 399 are 
replicated in 4U 2206H-54 - which is assumed to be powered 
by direct wind-fed accretion (Negueruela & Reig l200Tb - while 
the primaries in both systems resemble one another (Blay et al. 
[20061 Walborn et 

In each of the systems named above the exceptionally long 
pulsational periods of the neutron stars are explicable under the 
assumption that they had been significantly spun down, first via 
magnetic dipole radiation and subsequently the propellor mecha¬ 
nism (Li & van den Heuvel [19991 Reig et al. 1200911201 21 Popov 
& Turolla 120121l . Operating on comparatively short timescales 
(~ 10"^ - lO^yr; Popov & Turolla 120121 1. this would also be 
consistent with the relatively unevolved nature of VFTS 399 (cf. 
footnote 6 ), although it would require extremely high magnetic 
helds, > 10'"^G, in order to facilitate it (cf. Li & van den Heuvel 
|1999| and Appendix A). 

Nevertheless, under either scenario, it is of interest that 
VFTS 399 appears to host the second young pulsar within the 
30 Dor complex after the isolated neutron star PSR J0537-6910 
(~5arcmin/~ 70pc distant in projection). With Pspin ~ 16ms 
(Marshall et al. 119981 1. PSR J0537-6910 rotates much more 
rapidly than the putative pulsar within VFTS 399, naturally ex¬ 
plicable because of the lack of accretion driven spin-down in the 
former. However, if the VFTS 399 pulsar is exceptionally slowly 
rotating due to the effects of an extreme P-held one would be 
forced to explain their different field strengths {B > lO'^G ver¬ 
sus B ~ lO'^G for PSR J0537-6910; Marshall et al. [T998] l. de¬ 
spite the likelihood of both pulsars having formed from the same 
underlying stellar population (cf. Dufton et al. l201 111 . 

4.2. The 0-type star mass donor 

VFTS 399 demonstrates notable physical similarities with 
VFTS 102, another member of the high rotational-velocity O- 
type star cohort of 30 Dor (Ramirez-Agudelo et al. 120131 1. Both 
are emission-line stars, share the same spectral classihcation and 
their temperature, luminosity and surface gravity are identical 
within the errors quoted (Sect. 3.1 and Table 1 of Dufton et al. 
1201 lb . Based on the proximity of VFTS 102 to the isolated pul¬ 
sar PSR 10537-6910 (Marshall et al. rT998i) . Dufton et al. (|2^ 
explored a common binary evolution scenario for both objects, 
with the I 6 M 0 H- I 5 M 0 mass-transfer model of Cantiello et al. 
(12007b yielding a rapidly rotating secondary well matched to 
VFTS 102 at the point of supernova. 

While it is tempting to adopt a comparable pre-SN evolu¬ 
tionary scenario for VFTS 399 - albeit with the binary remaining 
bound at SN - this would be hard to reconcile with the extant RV 
data. Specihcally, a short orbital period would require the binary 
system to lose a large fraction of the mass and angular momen¬ 
tum released in the pre-SN mass-transfer process (Petrovic et 
al. lM)5ll . since conservative evolution leads to orbital periods 
in the range ~ 50 - 300 day (Wellstein et al. 2001). In other 
words, this observational constraint implies that VFTS399 ac¬ 
creted comparatively little material, noting that the accretion of 
only ~ 2 M 0 is sufficient to spin it up to critical rotation. Had the 
initial mass ratio been close to one - as inferred for VFTS 102 
- then the mass transfer would be expected to remain close to 
conservative (Langer l2012l de Mink et al. 12013b . Therefore the 


properties of the VFTS399 binary system are best understood 
if the initial mass of the NS progenitor was significantly higher 
than the current mass of the Oe component. 

Given the current absence of a full orbital solution for 
VFTS 399 the calculation of a tailored evolutionary model is 
premature, although based on previous analysis we might expect 
an initial mass of > 25Mq for the NS progenitor; as a conse¬ 
quence of this we may speculate that VFTS 399 might host a 
high-mass NS. Under the mass-transfer scenario we favour the 
mass gainer is expected to quickly rejuvenate, and then behave 
just as a single star of its new mass, with the possible exception 
of rapid rotation and surface abundance anomalies. Therefore, 
we might expect VFTS 399 to have been born with broadly sim¬ 
ilar properties to those we infer for it now (Table 5 and Sect. 3.1), 
while our current age estimate (~ 6 Myr; Table 5) would be only 
slightly lower than expected. Our hypothetical ~ 25Mq primary 
would be expected to explode after ~ 7Myr (Brott et al. 1201 lb : 
fully consistent with such a scenario and implying a rather recent 
(< IMyr) SN. 

As with VFTS 102, we might therefore expect the sys¬ 
temic velocity of VFTS 399 to have been influenced by the SN. 
Somewhat surprisingly, there is no indication of it being a ra¬ 
dial velocity runaway, with Viad ~ 289.4 + 18.6kms“' in com¬ 
parison to a cluster mean of ~ 271km ' and a velocity disper¬ 
sion cr ~ lOkms ' (Sana et al. 12012112013b . Nevertheless, the 
location of VFTS 399 - on the periphery of the complex and 
displaced from any recognised cluster (Fig. 1; Walborn et al. 
12014b - suggests a runaway nature; our prediction of a recent 
SN would suggest a displacement of < 50pc from its birth site. 
Ongoing Hubble Space Telescope proper motions studies (PI: 
D. J. Lennon; see Sabbi et al. 12013b should provide a comple¬ 
mentary determination of the transverse velocity in order to test 
this prediction. 

Both VFTS 102 and VFTS 399 have been associated with 
a population of apparently single, rapidly-rotating O-type star 
runaways within the 30 Dor complex (Walborn et al. 12014b . Of 
these, 15 stars are of spectral type 08-09.5 III-VPl and hence 
might be expected to exhibit the OeBe phenomenon (Negueruela 
et al. 12004b . However, at the time of our spectral observations, 
only VFTS 102 and 399 were found to support the gaseous cir- 
cumstellar discs from which a putative neutron star/black hole 
might accrete. Therefore the remaining stars from this grouping 
could also harbour relativistic companions which would remain 
undetectable at X-ray energies until they too formed circumstel- 
lar discs. Eldridge et al. (1201 lb examined the pre- and post-SN 
evolution of massive binaries and concluded that they remain 
bound in ~ 20 % of cases, with a most likely systemic velocity 
of ~ 50 kms ^ Given this prediction, long term optical and X- 
ray observations of this cohort would appear warranted in order 
to determine if a subset of these stars are indeed quiescent X-ray 
binaries and hence received their anomalous rotational and/or 
radial velocities via binary interaction and subsequent SN kick. 

5. Concluding remarks 

The multi-wavelength properties of the outlying, rapidly rotat¬ 
ing O-type star VFTS 399 appear difficult to reconcile with those 
expected for a single star, instead favouring a binary hypothesis. 
The X-ray properties strongly favour the presence of a relativis- 

“ VFTS 5 (08 Vn), 74 (09 Vn), 91 (09.5 Illn), 102 (09: Vnnne-t), 
138 (09 Vn), 249 (08 Vn), 399 (09 Illn), 530 (09.5III:nn), 
531(09.5III:nn), 574 (09.5 Illn), 592 (09.5 Vn), 654 (09 Vnn), 660 
(09.5 Vnn), 768 (08 Vn) and 843 (09.5IIIn). 
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tic accretor, with the optical and near-IR properties consistent 
with classihcation as a Be X-ray binary, being the seventeenth 
such system within the LMC at the time of writing (Negueruela 
& Coe 120021 Vasilopoulos et al. 120131) and the hrst identihed 
within 30 Dor. 

The discovery of an apparently robust, statistically- 
signihcant periodicity in the X-ray lightcurve provides strong 
evidence for a neutron star in the the system. Unfortunately 
we are currently unable to determine conclusively whether the 
~ 2567 s period reflects Pspin or instead is an alias of the true 
f^spin ~ 3.245 s, although the robustness of the former period 
to the secular evolution common in X-ray binaries argues that 
it is astrophysical in origin. Breaking this degeneracy via obser¬ 
vations with an alternative observatory such as XMM would be 
of considerable interest, since the former scenario would imply 
that the neutron star supported a magnetic held of comparable 
strength to those exhibited by magnetars in order to permit the 
dramatic spin down required to yield its current Pspin- 

These hndings have a number of ramihcations: 

(i) Identihcation of individual (high-mass) X-ray binaries with 
their natal stellar aggregates is surprisingly rare - LS I -1-63 235 
lies within the halo of the cluster NGC663, while PSR B1259- 
63 and IGR J00370-I-6122 are likely members of the Cen OBI 
and Cas OB8 associations respectively (Negueruela et al. l20111 
Gonzalez-Galan et al. 120141) . The location of VFTS 399 within 
the conhnes of the 30 Dor complex could therefore aid the ob¬ 
servational assessment of the predictions of current population- 
synthesis models. These include the production (and retention) 
rates of high mass X-ray binaries as a function of star-formation 
rate and the duration of such activity. This would also inform 
studies of e.g. the global X-ray luminosity/star-formation rate 
relationship for non-active galaxies (e.g. Mineo et al. 120121) . 

(ii) In conjunction with determinations of both radial and trans¬ 
verse velocity components, characterisation of the orbital period 
and eccentricity of VFTS 399 should permit the dynamics of the 
SN explosion to be reconstructed. Given that Pfahl et al. (120021 ) 
and Podsiadlowski et al. (120041 ) suggest that the magnitude of 
kick is dependent on the pre-SN binary evolution, the fact that 
we may place observational constraints on this for VFTS 399 via 
comparison to the stellar population(s) of 30 Dor is of consider¬ 
able interest. 

(iii) If the inference of an extreme P-held for the putative neu¬ 
tron star is conhrmed, reconstruction of the pre-SN evolution¬ 
ary pathway and SN dynamics will help constrain the formation 
requirements for such objects. Moreover, it would support the 
hypothesis that binarity is a common feature in the formation of 
magnetars, which in turn would arise from progenitors with a 
wide mass range from about 17Mo up to 4044© or more (Sect. 
4.2 and Clark et al. lma . 

(iv) After VFTS 102, VFTS 399 is the second object drawn from 
the population of rapidly rotating apparently single O-type stars 
for which historical/current binarity has subsequently been in¬ 
ferred. Such a hnding is consistent with the hypothesis that bi¬ 
nary interactions may play a key role in the formation of this co¬ 
hort (de Mink et al. 120131 Ramfrez-Agudelo et al. 120131) . If con¬ 
hrmed for further examples it will have considerable import for 
the production of the rapidly-rotating, chemically-homogeneous 
stars inferred to be the progenitors of long gamma-ray bursts. 
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(Reig et al. 120091 Sect. 2.2), which leads to an order-of- 
magnitude estimate of Z? ~ 10*"^G for the neutron star. One might 
assume that a signihcantly lower accretion rate may alleviate the 
requirement for such a B-held. However, following Li and van 
den Heuvel (119991 1 e.g. M — lO'^g s * (leading to B = lO'^G) 
results in a spin-down timescale comparable to the total lifetime 
of VFTS 399 (cf Dufton et al. 120111) . in contrast to its cuiTent 
relatively unevolved nature. 

Recently, Shakura et al. (120121 ) proposed a new theoretical 
model for quasi-spherical accession onto slowly rotating neu¬ 
tron stars (see also Popov & Turolla 120121) . Reig et al. (120121) 
demonstrate that the accretor in 4U2206-I-54 still supports an ex¬ 
treme (~ 10 '"^G) B-held under this formulation, although this 
requires an anomalously slow primary wind (which is observed 
in this system; Ribo et al. 120061) . Clearly, an orbital solution for 
VFTS 399 and a determination of the wind parameters of the 
primary (mass loss rate and velocity at the orbital radius) in or¬ 
der to determine M would provide an empirical test of the ex¬ 
treme B-held hypothesis, as would a measure of the long term 
spin-period history of the putative neutron star, from which an 
estimate of B may also be obtained. 


Appendix A: Spin-down in the presence of extreme 
magnetic fieids 

Following the analysis of Li & van den Heuvel (119991 ). the evo¬ 
lution of the spin period of the putative neutron star companion 
to VFTS 399 depends on the accretion rate/torque - a function of 
the (currently unknown) orbital separation and wind properties 
of the primary - and its magnetic held. For this mechanism to 
yield Bspin ~ 2567 s it must have been born with an (external) 
magnetic held of comparable strength to those exhibited by 
magnetars {B > 10^"^G). If we make the assumption that the 
neutron star in VFTS 399 is rotating at its equilibrium period we 
may estimate the dipolar held strength from; 

Bspin - 

(Li & van den Heuvel 119991 Reig et al. 120091) where B - 
10 ^^Bi 2 G is the dipolar held strength, M - lO'^Misg s ' the 
mass accretion rate and R = lO^Becm and M - L4Mi 4 M 0 the 
neutron star radius and mass. We adopt Pspin ~ 2567 s, R(, - 
Ml 4 = 1 and M 15 a; LxRjGM - 0.25 for ~ 0.5 x lO^^erg 
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